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Abstract 

We argue that, contrary to previous claims, the supersymmetric sine-Gordon model 
with boundary has a two-parameter family of boundary interactions which preserves 
both integrability and super symmetry. We also propose the corresponding boundary 
S matrix for the first supermultiplet of breathers. 



1 Introduction 



Ghoshal and Zamolodchikov IQ formulated the general framework of integrable quantum 
field theory in the presence of a boundary. Q Among the examples which they investigated 
is the boundary version of the sine-Gordon (SO) model |0, |[. They argued that this model 
has a two-parameter family of boundary interactions which preserves integrability, and they 
proposed the corresponding soliton boundary S matrix. 

Inami et al. subsequently investigated a boundary version of the supersymmetric sine- 
Gordon (SSG) model |0-|]12|. They argued that the combined constraints of integrability 
and supersymmetry do not allow any free parameters in the boundary interaction. The 
boundary SSG and SShG (supersymmetric sinh-Gordon) models have been studied further 
in ll-lpl. 



We argue here that, contrary to the claim in 0, the boundary SSG model has a two- 
parameter family of boundary interactions which preserves both integrability and supersym- 
metry. The key point is that the class of boundary interactions considered in is too 
restrictive: one must instead introduce a Fermionic boundary degree of freedom, as was 
done by Ghoshal and Zamolodchikov fl]] in order to describe the Ising model in a boundary 
magnetic field. The boundary SG model at the free Fermion point [l7, 18] is also formulated 
with such Fermionic boundary degrees of freedom. 



Our result resolves an interesting paradox which arose from our recent work on the 
boundary supersymmetric Yang-Lee model, which is defined using the notion of perturbed 
conformal field theory PU|. Indeed, let us first recall that, just as the bulk Yang-Lee (YL) 



model 1^ is a restriction p2| of the bulk SG model, the bulk supersymmetric Yang-Lee 
(SYL) model ||2^ is a restriction of the bulk SSG model. The same is true for the 
corresponding boundary versions: the boundary YL [|l], ^ and boundary SYL |jT9| models 



are restrictions of the boundary SG and boundary SSG models, respectively. In we 
found strong evidence (both from perturbed conformal field theory and the conjectured 
boundary S matrix) that the boundary SYL model with one free boundary parameter has 
supersymmetry. This strongly suggests that its "parent", the boundary SSG model, should 
also have at least a one-parameter family of boundary interactions which maintains both 
integrability and supersymmetry. 

The outline of this Letter is as follows. In Section 2, we propose the boundary interaction 
in terms of two unknown functions /(0) and B{(j)). We determine these functions in Section 
3 from the requirements of supersymmetry and integrability. In Section 4 we briefly discuss 
our results. In particular, we propose the boundary S matrix for the first supermultiplet of 
breathers. 



^Earlier related work includes (|]-|^. 
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2 The boundary SSG model 



In order to facilitate comparison with 0, which hereafter we refer to as I, we adopt the same 
notations. The Euchdean-space action of the boundary SSG model is given by 

/oo /"O /»oo 

dy dx Co + dy Cb , (1) 
oo J — oo J — oo 

where the bulk Lagrangian density is given by^ 

£o = 2d,(j)d2(f) - 2^d^^ + 2^8^^ - 4 cos - A^i) cos ^ , (2) 

where if) and are the two components of a Majorana Fermion field, and z = x + iy, 
z = X — iy. We propose the following boundary Lagrangian at a: = 

Cb = iJ^ + iadya - 2f{(f))a{^ - ^) + B{4>) , (3) 

where a is a Hermitian Fermionic boundary degree of freedom which anticommutes with 
both ifj and if). As mentioned in the Introduction, such an approach was first considered 
by Ghoshal and Zamolodchikov to describe the Ising model in a boundary magnetic 
field. Moreover, /(0) and are potentials (functions of the scalar field 0, but not of its 
derivatives), which are still to be determined. 

Variation of the action gives the classical equations of motion. For the bulk, the equations 
are (12.2), (13.3) 

^5(9^0 = sin + ^^^Z; sin ^ , 

dzip = -i> cos ^ , d^ip = -iIj cos ^ , (4) 
and the boundary conditions at a; = are 

^ + ^ = 2/a, (5) 
idya = f{ip - ip) , (6) 

If we assume that / is nonzero, then we can eliminate a from Eqs. (^-(0)) and thereby 
obtain the following boundary conditions at x = 0: 

^y{^P + ^lJ) = ^dy<P{ij + ij)~2tf{ij-i;), (8) 
r. , dB d\n f - , 

d^(j) = - +2—ftljil) . (9) 



^We have already performed, following I, the field rescalings (t> ^ <t>IP, ip ^ ip ~^ V'//?, and we 

consider the classical limit /3 ^ 0; moreover, we have set the mass m = 2 (cf. Eq. (12.1)). 
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3 Integrals of motion 



We now proceed to determine the potentials and B{(f)) which appear in the boundary 
action by demanding both supersymmetry and integrabihty. We recall |11|] that the bulk 
SSG model has an infinite number of integrals of motion constructed from the densities Tg+i, 
T^+i, 9^_i, Qs-i, which obey 

d,T,+i = d,Qs-i , dJTs+i = dM~i ■ (10) 
The s = I and s = 1 densities are given by the supercurrent 

e_i = -2^sin^, G_i=-2^sin^, (11) 

and the energy-momentum tensor 

T2 = {d^cPf - d^ij^ , % = {d,(Pf + dd^ , 

©0 = 00 = — 2 COS — V'-i/' COS ^ , (12) 
respectively. The s = 3 densities are given by P, pT]] 

T, = {^l<P)'-^{^.<py + ^{^.<P)'^.^P^^;-^l^P^.^, 

T, = (a|0)2-^(a,<^)^-^(a,0)2a,^^ + 52V'a,V^, 

©2 = {dz(py COS (p - ^z^p^^J cos^ ^ - ^^J^z^|J^;, COS ^ + ^^ljijj{dz(j))'^ COS ^ , 

02 = {^,(j)ycos(j) + ^,iP^Pcos^^ + iJ^,^P^,cos^ + ^iPiJ{^,(PYcos^. (13) 

As observed by Ghoshal and Zamolodchikov , it follows from the continuity Eqs. (|T0|) 
that the boundary model has the integral of motion 

Ps= [ dx {Ts+i + T,+i + e,_i + 0,_i) - iT.,{y) , (14) 

provided that the following condition holds at x = 

T,+i - T,+i - e,_i + e,_i = dyj:M , (15) 

where Ss(y) is a (local) boundary term. 

As we shall show below, the requirement that P3 be an integral of motion implies that 
the potential B{(t>) is the same as for the (Bosonic) boundary SG model 

B(0) =2acos^(0-0o), (16) 



where a and 0o are arbitrary real parameters. 

The requirement of on-shell supersymmetry is that Pi also be an integral of motion. The 
constraint (^) for s = ^ reads 

{}^dA - 2 sin ^){^ -^)- '-dy(t>{i^ + V^) = 1 . (17) 

We now substitute for d^ct) the boundary condition with the potential B{(j)) given by (]T6|); 
and then we use (|^) and (||) to substitute for ■?/' + ■0 and i/j — ip, respectively. Moreover, we 
make the following Ansatz for the boundary term 



^ = -g{<P)a, (18) 

where a is the Fermionic boundary degree of freedom appearing in the boundary Lagrangian 
(^), and g{(f)) is a function of which is yet to be determined. Eq. ( P7| ) implies that g{(f)) 
must satisfy the pair of equations 



a sm 



— 4 sin 



dg 



f 



-2/. 



(19) 



It is not difficult to solve these equations for /(0) (and hence, also g{(p)), 



a sin - ( 



(po) - 4 sin - 



a cos ■ 



0o) -4cosf + C 



(20) 



where C is an integration constant. 

We now turn to the requirement that P3 be an integral of motion. The LHS of the 
constraint (^) for s = 3 is given by (13.10) 

r4 - T4 - 62 + 62 = 



-(4 sin (f) — 2(f)yy + 2ifjilj sin 



+ 



hy + i(f)x<Py COS ( 



3i . . 

+ - (fi)i'^y^ - ^yi^) + 3« COS^ ^ii^yip - ^yip) + ^ sin ^MipyiJ - ipytp) 

(j) 3 

+ i{'ljjyyll!y " i^yyi^y) + ii)ll)(j)x(t)y COS - " - (f) ^ (f) y {l/j ylj) + l/jyl/j) 

Z o 



+ 2dy COS -{ipyi^ - Ipi^y) + COS -{ipi^yy " i^yyi^) ■ 

Our task is to try to write the above expression as dyT,-^. 

We begin by replacing (j)^ by the boundary condition (^, and (j)xy by 



^xy 



(21) 



(22) 



Since the pure Bosonic terms in (^) (which evidently do not depend on the potential /) 
must by themselves form a total ^/-derivative, it is clear that B{(f)) must be the same as for 
the pure Bosonic boundary SG model, as anticipated in Eq. (|1^) above. (See also I.) 

The last two terms in (^) can be re-expressed as 



3dy cos ^{ipyip - V'^y) + dy 



cos -{ipl/jy - l/jy^) 



(23) 



Let us now consider the Fermionic term in (^Tj) with three derivatives, i{iljyyifjy — ipyyipy) 
We can extract a total y-derivative from this term as follows: first, observe the identity 



dy{lpy^y) = i^yy^y + i^y^^yy ■ 

Next, differentiate the boundary condition (|]), 



(24) 



"^yy ~^ "^yy 
dlnf 



k2 ^y 



+ ^rr^yy (V' + ^) 



d 111 f 



(25) 



Solving 
result 



]) for i/jyy and substituting into (p4|), and repeating for ipyy, we obtain the desired 



ii^l^yy'^y " ^yy^y) = -idy^^Jy^y 

'a^ln/. 51n/, 



k2 ^y 



+ 



^yy 



^y 



+ 2l (PyWyWy ■ 



(26) 



Our general strategy now is to use the boundary condition (|^) to express ifjyilj in terms of 
ipyif} (plus terms with ipip)] and similarly to express ifjyip in terms of ipyip-, and ipy'ipy in terms of 
ipyif} and ipyip- In this way, all the remaining Fermion bilinears have at most one y-derivative; 
and those bilinears with one derivative appear in only two combinations: (^/'y'0 — 4'y'4^) and 
('?/'y'0 -|- ipyifj). The former combination can be further reduced using the identity 



dlnf - 

Ipyip - Ipyl/J = dy^lpl/j) " 2 (j)y1p'l{) 



(27) 



However, the latter combination cannot be so reduced, and thus, its coefficient must vanish. 
After some computation, we find that the total contribution of such terms in (p^) is 



{iPytp + 1pytp)^y 



a sin - ( 



bo) - 4 sin - 



+ 12/ 



,d\nf 



{2i 
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Recalling the result pO|) for /(0), we find that the expression in (p8[) within braces does 
vanish, provided that the integration constant C is given by | 



C 



o? — 8a cos h 16 



(29) 



The expression for then takes the simplified form 



C . 1 
— sm - ( 
2 4 



where tan 



D 



a sm ^ 



2 a cos ^ — 4 



(30) 



After further computation, we find that (^) is given by 

3 . , 3 . (/Jn 3 



- sm - -a sm y - —a sin(0 - ^o) - 6/ sin - 



+ I ^a^sin^ ^1 



6o) + 6 cos - 6 cos^ I - 24/ 



2\ 51n/ 



a^ln/ ain/a^ln/ 1 <91n/ (dXxvj 



16 



3 aMn/ /<91n/ 
16+^^^ + ^ 



(31) 



up to a total y-derivative. Remarkably, upon substituting the result (|30|) for /(0), we find 
that the quantities in each of the square brackets vanish. We conclude that P3 is indeed an 
integral of motion. 

Finally, we remark that the boundary model also has the integral of motion 



K = [ (r,+i-T,+i + e,_i-e,_i)-2s;(i/) 



provided that the constraint 

Ts+i + Ts+i - Qs-i - Qs-i = dy^'s{y) 



(32) 



(33) 



is satisfied at x = 0. We expect that the boundary SSG model has (for s = |, at least) such 
integrals of motion if a different set of boundary conditions is imposed, corresponding to a 
boundary Lagrangian of the form 



= -ipip + iadya - 2z/(0)a(V' + ip) + B{(p) . 

However, we shall not pursue those calculations here. 

•^There is a second (negative) root, which is ruled out by the requirement that / be real. 



(34) 
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4 Discussion 



We have shown that the boundary SSG model with boundary Lagrangian (|^), with B{(j)) 
given by (|T6|) and given by (|30|), has the integrals of motion Pg (|T3|) with s = | and s = 3. 
The fact that Pi is an integral of motion means that the model has on-shell supersymmetry. 
The fact that P3 is an integral of motion is strong evidence that the model is integrable. 
We emphasize that this boundary interaction, just like the one for the (Bosonic) boundary 
SSG model, has two continuous parameters, a and 00 ■ As mentioned in the Introduction, 



we were motivated to look for such boundary interactions by our earlier results |jT9[ on the 
boundary SYL model. 

We conjecture that the boundary SSG model is integrable, and that the corresponding 
boundary S matrix S(^) for the n = 1 breather supermultiplet is given by |jl5| 

S(^) = SsGiO;v, ^) S^susY iO;ip), (35) 



where Ssg{(^ '^V ^'^) is the n = 1 breather boundary SG reflection factor p6[, and S^susvi^ j V^) 
is a 2 X 2 matrix given in [^, |T3|. Indeed, S(^) satisfies the boundary Yang-Baxter equation, 
as well as unitarity and boundary cross-unitarity Moreover, S(^) commutes with a 
supersymmetry charge corresponding to the integral of motion Pi .0 The proof is essentially 
the same as the one given in for the boundary SYL model. Note that the boundary term 
Si ([T8|) corresponds (up to a proportionality factor) to (—1)^, where F denotes Fermion 
number. 

Although the 5* matrix (^) involves three boundary parameters r/, we expect that 
there is one relation among these parameters, so that only two parameters are independent, in 
concordance with the boundary action. Indeed, a parallel situation occurs for the boundary 



SYL model [0, whose boundary S matrix nominally involves two boundary parameters, 
but only one of them is independent. For the boundary SSG model, the relation among the 
parameters can presumably be found in the same manner as for the boundary SYL model 
]TP| : namely, by imposing the constraint []T| that near a pole iv^^ of the boundary S matrix 
associated with the excited boundary state la)^ (which can be interpreted as a boundary 
bound state of particle Aa with the boundary ground state |0)b), the boundary S matrix 
must have the form 

S^(g) - , (36) 

where g'^Q are boundary-particle couplings. It would be interesting to explicitly work out this 
relation for the boundary SSG model, and to find the precise relation of the two independent 
parameters of the boundary S matrix to the parameters a, of the boundary action. The 
latter problem has already been addressed for the case of the (Bosonic) boundary SG model 



27| , |28[| , and we expect that similar techniques can be applied to the supersymmetric case. 
Moreover, analogous results should also hold for boundary S matrices for the SSG higher 
breathers and solitons. 



"^The matrix S^^^^y(0 ; Lp) with e = — 1 commutes with a supersymmetry charge corresponding to Pi. ( |3^ ) 
.1, „ _ 1 



with s = \ 
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Note added: 



We have verified that the boundary SSG model with the boundary Lagrangian CJ^ ( p4D 

indeed has the integrals of motion P[ (|32D and P3 (p!^ , with the same (|^), and 

2 

= ^cos|(0 — D) where C and D are given by (p9D and (0) with a —* —a. In the 
conformal limit, the boundary Lagrangians and (^4|) give Neveu-Schwarz and Ramond 
boundary conditions, respectively. 
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